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ABSTRACT 



The problem of determining whether a line of bearing 
measured by a local surface-based sensor coincides with a 
threat whose position has been previously estimated by an 
ocean surveillance sensor is addressed. Uncertainties in 
the position estimate of the threat, in the bearing error 
and in the position estimate of the sensor are considered 
in measuring the probability that the threat lies on a given 
bearing from the sensor. A TI-59 calculator program is 
developed which calculates this likelihood when the threat 
location density can be assumed to be bivariate normal. 
Computations required when significant time has elapsed 
since the original estimate of threat location when the den- 
sity can no longer be considered bivariate normal are 
discussed . 
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INTRODUCTION 



I . 

As long range surface launched weapons systems continue 
to be introduced into the fleet, the operational commander 
is increasingly faced with the problem of being able to 
launch weapons at targets located beyond the horizon. 
Successful employment of such weapons is dependent not only 
on the ability to detect, classify and localize targets at 
considerable distances, but also on the ability to distinguish 
the true target from a potentially larger field of false 
targets. While long range ocean surveillance sensors may be 
of assistance in the identification and localization oi tar- 
gets, the information provided may not be refined sufficiently 
to permit effective targeting of long range weapons on that 
basis alone. The on-scene commander must in general rely on 
additional data on target location gathered locally and close 
to the time of weapons launch for accurate targeting. Thus, 
he must still be able to detect and track the desired target 
and be able to distinguish it from other targets within 
range of his sensors. 

The procedures developed in this paper are designed to be 
of assistance in addressing the last of these problems. They 
are applicable when the information available is an error 
ellipse around a threat location estimated by an ocean sur- 
veillance sensor and bearings only data generated by a local 
surface-based sensor. The question of whether a target 
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detected by the local sensor is the same as that whose 
estimated position was provided by an external sensor can 
only be addressed if information is available on the loca- 
tions and tracks of all possible targets within range of the 
local sensor. Since such data is generally not available, 
this paper does not attempt to answer that question, but 
rather develops a method by which bearing information from 
different sensors can be compared as to the likelihood of 
each bearing being associated with the threat identified 
previously. It is envisioned that these likelihoods can 
be then used to induce an ordering among bearing data gathered 
by different sensors or, conceivably, conflicting data 
gathered by one sensor. The ordering would be based on the 
likelihood that each bearing will contribute to refining the 
original estimate of the location of the target of interest. 
This information could be applied in a number of ways: as 

a guide to allocation of more capable sensor resources for 
purposes of obtaining targeting information; as a guide for 
allocation of weapons against more than one threat; as a 
means of pre-processing data before entering it into a target 
motion model, thereby reducing the chance of introducing 
unrelated data. 

To determine the likelihood that a given line of bearing 
and the threat coincide, consideration was given to the 
uncertainties inherent in estimation of target position, 
in the measurement of bearings by a particular sensor and 
in estimation of sensor location. It is assumed that at some 
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time t , an ocean surveillance sensor detects a threat whose 
o 

position is estimated to be within an elliptical region with 
p^ x 100% certainty. The estimated position data are received 
and converted by the on-scene commander into a probability 
distribution described by a truncated bivariate normal density 
function. 

It is further assumed that the standard error a. charac- 

3 

teristic of the local sensor is known. The sensor bearing 3 
with bearing error a. is then projected from the sensor 
position through the threat density function. 

Since sensor position relative to the target may itself 

be subject to navigation error, the uncertainty is introduced 

as a truncated circular bivariate normal distribution centered 

at location (u , v ) with standard deviation a. 
o o 

A TI-59 calculator program is developed which estimates 
the likelihood that the threat identified by an external 
sensor lies along bearing 3 , given the threat distribution, 
the bearing error, and sensor position distribution relative 
to the threat. 

The theoretical basis for this calculation is presented 
in Chapter II. The algorithms used in designing the calcu- 
lator program are described in Chapter III. A program listing 
and verbal flow are provided in Appendix A along with instruc- 
tions for the user. Appendix 3 contains a verbal flow of a 
program designed for use when considerable time has elapsed 
since the initial estimate of the location of a moving 
target . 
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II . THEORETICAL BASES 



The general approach to determining the likelihood that 

a measured line of bearing 3 is the true bearing from the 

sensor to the threat identified and localized by an external 

sensor is discussed in this chapter. Calculations required 

when using threat position information both as initially 

generated by the external sensor at time zero, t Q , and as 

distorted to account for an intervening time late, t , are 

J-j 

discussed. Uncertainties in bearing measurement and sensor 
position are included. 

Initially, assume that sensor position is known with 
certainty. Let 3 T be the true bearing of the threat from the 
sensor. Since threat location is uncertain, 3 T is a random 

quantity with probability density function f (3 ) . Let 

; TRUE 

f (3;3 rp )d3 be the probability that the errors in bearing 
measurement are such as to give rise to a bearing on the 
threat in the interval dB about the observed value 3 when 
the true bearing is 3 . Then the likelihood of observing a 
bearing 3 is: 

f'(6) = / f Q (B; Bm) f R (6m) dBm • (D 

all B t * T "TRUE 1 

The density f (3 ) is determined by the probability 

~TRUE T 

density function of threat location at time t, f v (x,y;t). 

A , X 
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Given the probability density function of threat location 
at time t and assuming that bearing errors are normally dis- 
tributed with mean zero, consideration is limited to the 
probability that the observed bearing is the true bearing 
of the threat, given the threat is contained in a planar 
region E and the true bearing $ T is in an interval about 3 
with upper and lower bounds 3 and 3 . The region E is 
selected to be the minimum area planar region which contains 
the threat with a specified high probability p^ . The inter- 
val (3 T ,8 TT ) is selected so that, for all the possible values 
of 3 t contained in the interval, 3 is within an interval of 
specified high probability p, around 3 T - The fan 
is symmetric about 3 so that it represents the minimum area 
region which meets the above criterion. The likelihood of 
observing 8 when the threat is in E and S T is in (8^,8^ is 
determined by integrating expression (1) above, over a region 
defined in the manner of the shaded area of Figure 1 (a). 

If sensor position is not certain, assume that it is 
distributed in accordance with a circular bivariate normal 
distribution f^. ^(u,v). Again, consideration is limited to 
determining the likelihood of observing a bearing 8 given 
that the threat is contained in a planar region, E, 6m is in 
( ^ L ' ^ u ) ' an ^ t ^ ie sensor contained in a planar region C. 

The region C is selected as the minimum area planar region 
which contains the sensor with a high probability p~,. The 
region of E over which f 1 (3) is evaluated expands as the 
shaded region of Figure 1 (b). 
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(a) Sensor Position Certain 




FIGURE 1. REGION OF INTEGRATION 
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Since the measurement errors involved in estimating the 
threat position, the sensor position and the bearing angle 
arise from different measurement procedures, it is reason- 
able to assume that errors are independent. The probability 
densities required to perform the above calculations can be 
estimated as follows. 

At time t an ocean surveillance sensor estimates the 
o 

position of the threat to be located within an elliptical 
area characterized by the parameter set E = {X,Y,9,A,B} with 
confidence p, x 100%. The elements of the parameter set E 
are respectively: X the latitude of the estimated threat 

position, Y the longitude of estimated threat position, 

9 the orientation of the major axis from true North, A the 
length of the semi-major axis, and B the length of the semi- 
minor axis. Since the measurement errors in determining 
threat position are generally assumed to be normally distributed, 
the ellipse characterized by E represents the minimal area 
p^ x 100% confidence region about the mean (X,Y) . Treating 
this ellipse as the p^ probability region of a bivariate 
normal distribution, a density function for the threat posi- 
tion at time t can be estimated. For convenience, locate 
o 

the origin of a rectangular coordinate system at the center 
of the ellipse, (X,Y) , with positive x-axis located along 
the major-axis of the ellipse at a bearing 9 from true North. 
Assume a flat earth in the region of interest. Let t = 0. 

The mean of the threat position density f (x,y;0) is then 
the point (0,0). The variances in the X and Y directions can 
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be derived from the fact that the region with minimal area 
which contains the threat with probability p, is a k-sigma 
ellipse where k is determined from the relationship: 

_ k 2 

P [threat located in k-sigma ellipse] = 1 - e 



Thus , 



a x 2 = (A/k) 2 

and 

a Y 2 = (B/k) 2 



where 



k = »-2 ln(l-p^) . 

X and Y are assumed to be independent. 

If the course and speed of the threat are known with 
certainty to be and s respectively, the probability density 
of the threat position at time late t. r can be shown to be 
again a bivariate normal with mean (st T cos ( 6-\p ) , st sin ( 9-^) ) 

1j Lj 

. 2 2 

and variances a , a . The p. probability region of the 

X ± -L 

density at time t would then be an ellipse congruent to that 
characterized by the set E above but centered at the point 
(st T cos ( e-ii/) , st T sin (e-i^) ) (Figure 2) [Ref. 1] . 

Li 
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FIGURE 2. TIME LATE ELLIPSE, KNOWN COURSE AND SPEED 
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In some cases, the motion of a submarine on patrol in a 

large area can be characterized, when t T is large, by an 

-Li 

expansion of the probability area with time at a rate D. 

The result in such a case is that f v v (x,y;t T ) is still 

A f x 

. 2 
bivariate normal with mean (0,0) but with variances a v + Dt 

A Xj 

2 

and a + Dt . 

Y 

When the motion of the threat cannot be described in the 

above manner, and the course and speed are not known, but 

assumed to be distributed according to the densities f (p) 

and f c (s), determining the probability density f^ Y (x,y;t„) 

is a considerably more complex problem. Let (x q (s , p) ,y (s, p) ) 

be the coordinates of the point at which the threat would 

have to be located at time zero in order to reach the point 

(x_,y T ) at time t T if the threat speed were s and course p. 

Thus, x (s,p) = x_ -st T cos (9-tjO and y (s,40 = y_ -st_ sin ( d-p) 

O L L Q Li Li 

(Figure 3) . Then, 



- 360 

f X/Y (x,y;t L ) = J J f XjY [x 0 (s,.»,y 0 (s,*);t 0 ] 

• f, F (p) f s (s) dtds [Ref. 1] . 

This density is no longer normal. In the special case when 
2 2 

A = B, i.e., a v = a , p has a uniform distribution over 
the interval (0°,360°) and S is known with certainty, as 
shewn in Reference 2, the density changes with time as in 
Figure 4. 
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FIGURE 3. TIME LATE POSITION, COURSE AND SPEED UNCERTAIN 
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FIGURE 4. TIME LATE THREAT DENSITY , SPEED KNOWN 
COURSE UNIFORM (0° / 360°) 
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The distribution of the bearing error measured by the 
local sensor can be estimated if the standard error of the 
sensor, a , is known. The bearing error is then assumed to 

p 

2 

be normally distributed with mean zero and variance c . 

3 

Assume that the estimated sensor location (u ,v ) is 

o o 

accurate to within R nautical miles with p^ x 100% confidence. 

The density f (u,v) of sensor location can be assumed to be 

U / V 

a circular bivariate normal with mean (u ,v ) and variance 

o o 

2 2 2 

a = R /k . The value of k is determined from the relation- 
ship : 

2 

_k /2 

P[(u,v) contained in k-sigma region] = 1 - e ' , 

where the probability on the left is p^ in this case. 

Since the evaluation of f'(8) considering sensor position 
density and bearing error density does not generalize to a 
closed form, algorithms are developed in the remainder of this 
paper for estimation of the probability (likelihood) that the threat 
identified by an external sensor lies on a line of bearing 
measured by a local sensor given that the threat is in region 
E, the sensor is in region C and the true line of bearing 
lies within 
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III. ALGORITHMS 



As indicated in the previous chapter, the variety of 
geometrical situations which can arise depending on the 
location of the sensor relative to the estimated threat loca- 
tion precluded development of a generalized analytical pro- 
cedure. Rather, algorithms are developed in this paper for 
numerically evaluating f'(3). The procedures applicable to 
a bivariate normal threat location density have been imple- 
mented on a TI-59 calculator. Appendix A contains a listing 
of that program. The calculations required when the time 
late threat location density is no longer normal exceeded 
the available program capacity of the TI-59 and therefore 
have not been implemented. A detailed verbal flow is provided 
at Appendix B for future implementation on a larger machine. 
The algorithms used in both situations are described in this 
chapter . 

A. BIVARIATE NORMAL THREAT LOCATION DENSITY 

This case includes situations (1) where the time elapsed 
since generation of the initial error ellipse by the ocean 
surveillance sensor is negligible, (2) where the motion of 
the threat can be assumed to be random in the manner described 
above, and (3) where course and speed of the threat are 
assumed to be known with certainty. With appropriate modi- 
fications to the input data, all three of these situations 
can be addressed using the program contained in Appendix A. 
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In situation (i) the data entered are the parameters of the 
ellipse as generated by the ocean surveillance sensor. In 
(2) , the location and orientation from North of the ellipse 
is the same as originally generated, but the size of the 
ellipse expands at some constant rate of area per unit time 
which must be estimated by the user. This rate D, times the 
elapsed time, t T , yields the factor which must be added to 
the semi-ma jor and semi-minor axes of the original ellipse. 
That is, if the original error ellipse is a p x 100% confi- 
dence ellipse, the semi-major and semi-minor axes of the 
diffused ellipse will be input as A' and B' respectively: 



In situation (3) the dimensions and orientation of the time 



late ellipse are the same as those of the original error 
ellipse, but the center of the ellipse is displaced from 
its original position by the known velocity times elapsed 
time. The updated position of the error ellipse is treated 
as the origin of the rectangular coordinate system for this 
situation and all linear measurements are made relative to 
this system. All angular measurements are made from true 
North. 

Estimation of the likelihood that the threat lies on a 
bearing from the local sensor given that the threat is located 



A' 




B' 



VB“ + (-2 In (1-p, ) ) Dt 

L L 
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in a p 1 x 100% confidence ellipse and the true bearing lies 
within the bounds (6 ,S TT ) with confidence p 0 x 100% proceeds 

L U Z 

as follows: 

1. Estimate parameters of the bivariate normal density 

2 2 

f X,Y^ X ' y ) of threat location: u x = u y = 0, <? x , a y . 

2. Determine 3 T and B TT such that an interval of length 

Li U 

2k a centered on either 3 T or 3 TT would contain 3 , the 

P P L U 

measured bearing, with probability p 0 : 3 T = 3 - k: _ a 0 , 

Z Li p p 

3 u =3 + Va‘ 

3. Determine sensor location coordinates relative to 
the origin of the threat ellipse. 

4. Subdivide the angular interval (3,3) into 2n sub- 

L U 

intervals. Each subinterval k intersects the ellipse in a 
strip with average width which corresponds to A3 fc . 

5. At the midpoint of each subinterval k determine the 
equation of the line through the sensor position at the true 
bearing 3 ^ from North. 

Y-c 

6. Let the equation of the line of bearing 3^ be X = 

Then the plane perpendicular to the xv-plane which contains 
this line intersects the bivariate normal threat density in 
a curve whose equation is 



g(y) 



1 

2 TTa x c x 




e 



(y-c) 2 



2 

m a x 



2 




( 2 ) 



Y-c 

found by making the substitution X = m the density 

function f v (x,y) . It will prove convenient to expand the 

A , I 

right side of (2) as follows [Ref. 3]: 
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1 

2 



2 



c 



g(y) 



m 



e 



r 2 2 2 
°Y m °X 



(3) 



o — / 2 2 2 

/2 :r \/a + m a 

Y X 




2 



/ 





e 



v 2 7T a v s v m 

X Y 



For computational purposes, assume that the width W v 
of the region of the ellipse cut out by the angular subinter- 
val around 3^ will be nearly constant through the ellipse. 

Let the points of intersection of the line of bearing 
with the threat ellipse be (X]^,Y^) and (X ,Y . Then 
approximate the volume of the normal density over this region 
by the absolute value of the product of the area under the 
curve g between Y^ and Y and W (Figure 5) . Observe that 
the term in brackets in equation (3) above is the density 
function of a univariate normal random variable with mean 



The term preceding the brackets in (3) is the slope of the 
line of bearing, m, times the density of a univariate normal 




and variance 




2 2 

random variable with mean zero and variance Cy + m 
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FIGURE 5 



evaluated at c. Thus, evaluating the expression 



|w k / g (y) dy| 

Y t , 



is equivalent to the computation 




1 



c 



;) ) (*(Z 2 ) - $(Z X ) ) | 



where 



Z . 

l 




for i = 1,2, i s the ~P[Z<_z] when Z is a standard normal 

random variable, and <p is the density function of a standard 
normal random variable. 

If the substitution Y = mX + c were made in the density 
f v v (x,y), an analogous situation would arise with the limits 

A , x 



of the integration being X,,, and X^. ~ . 

7. If the value of the density under the bivariate 
normal curve over the region of the ellipse defined by the 
subinterval AB^ is weighted by the probability of observing 
a bearing error (3^-3) the result is the probability of 
observing the bearing of the threat as 3 when the true threat 
location is in the segment of the ellipse defined by AB^- 
This probability is approximated as follows: 
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k2 

\ y I f x,Y ( 5 T'y )d yi V 

k1 



Since the bearing error is assumed normally distributed with 

2 

mean zero and variance a , the value of f (3; S,) can be 

P 3 K 

determined by the expression 



_L ) , 

a „ a „ 



where <p is as above the density function of a standard 
normal random variable. 

8 . 



f ’ (3) = f f (3 ; S t ) f 



all S 



TRUE 



<3 t > da t 



is approximated by 



2n 



k2 



I ! w >t / f x,Y (Z ^T' y)dy] I f ^ ( 3 ; 3 v) 



k=l 



kl 



The value of the sum is determined by repeating steps 
A. 6 and A. 7 above at the midpoints of each of the 2n sub- 
intervals defined in step A. 4 and summing the results of 
each of these calculations. Obviously, the finer the sub- 
division of (3 t /S tt ), the more accurate will be the estimate 

1j u 

of the likelihood, but also the longer the calculation will 
take. 
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B. SENSOR POSITION UNCERTAINTY 



The result of the above calculation will be the likeli- 
hood f ’ ( 3 ) that the threat lies along the bearing measured 
by the sensor given that the threat is located within the 
threat ellipse and the true bearing is within the interval 
(S^/Sy) and the sensor is at the position used to perform 
the calculation. We next will introduce additional calcu- 
lations that are required to account for the fact that the 
sensor position is not known with certainty. 

1. If the position of the sensor is estimated as being 

within R nautical miles of ( u 0 / v 0 )' its assumed coordinates 

in the xy-system, with p , x 100% confidence, estimate the 

parameters of the sensor location density f ^(u,v) with 

2 

mean zero and variance a : mean = (u , v Q ) , variance, 

2 2 

a = R / ( — 2 ln(l-p^)) in threat centered coordinates. 

2. The bearing measured by the sensor is B regardless 

of the sensor location. Assume that the area of intersection 
of the angular wedge (3^,3^) and the threat ellipse does not 
change significantly as the sensor position is moved along 
the line of bearing 3 . Then the effect of the bivariate 
normal distribution of sensor location can be approximated 
by considering only the univariate normal density along a 
line through (u , v ) perpendicular to 3. Repeat the calcu- 
lations in steps A. 3 through A. 8 above with the sensor located 
at each of the three points (u ,v q ) , (u q + . 97gcos ( 0-3-90) , 

v +. 97asin ( 6-3-90) ) , and (u -. 97ocos (9-3-90) ,v -97asin ( 8-3-90) ) . 
o o o 

If the line through (u ,v Q ) perpendicular to 3 is subdivided 
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symmetrically about (u , v ) such that 1/3 of the univariate 

o o 

normal density lies above each subinterval, the three points 
chosen above represent the "center of gravity" of each third 
of the density (Figure 6) . 

3. If R is chosen to include a significant proportion 
of the sensor density, i.e., on the order of 2 a or greater, 
the probability of the sensor being located in each of the 
three regions is approximately 1/3. Thus, if p, is on the 
order of .86, multiply each result in step B.2 by 1/3. 

4. Summing the results of steps B.2 and B.3 yields an 
estimate of the likelihood that the threat lies at bearing 
8 given the threat is in the p^ * 100% confidence ellipse, 

3 t is in and the sensor is in a p... x 100% confidence 

region. That is, the likelihood is estimated by 

3 1 

f ' (3) = S 4 f . ' (3) , 

j-1 3 

where j is the index of sensor position in figures 5 and 6. 

Instructions for application of the TI-59 program to 
calculate the above are included in Appendix A. 

C. THREAT DISTRIBUTION NOT BIVARIATE NORMAL AFTER TIME 
LATE ELAPSED 

The basic approach to evaluating f'(3) when the time late 
distribution of the threat is no longer bivariate normal 
is the same as that just discussed. The principal difference 



28 



N 




v o -.97csin( 0-3-90) ) v q +.97 sin ( 0- ^" 90 ^ 

= ^ u 3 ,v 3^ ~ ^ u 2' v 2^ 



FIGURE 6. ESTIMATE OF SENSOR POSITION DENSITY 



29 



arises from the fact that the time late density is signi- 
ficantly more complex in this case. 

The threat distribution becomes distorted from the normal 
after some time late when the course and speed of the threat 
are constant but not known with certainty. Application of the 
method described herein requires that the user assume a dis- 
crete distribution of the speed of the threat with upper and 
lower bounds S,^ and S^, respectively. In addition, the threat 
course is assumed to be uniformly between 0° and 360°. 

The density of threat location after some time late t_ 

Li 

when the threat speed is s^ then becomes [Ref. 1] : 



f X,Y <X,y;t L' S i ) = 



, 360 , x-s . t cos ( ) 

1 ( r J. / 1 Li 

=5 exp [ -7T ( - — - — - 

27Ta,, a,, c 2 2 



X Y 0 



X 



(y-s.t L sxn(e-*)) ^ diJ 



360 



, j 

where 9 is the bearing of the major axis from North. Note 

that the new threat density is still centered at the same 

position as the time zero ellipse but its shape changes as 

in Figure 4 of Chapter II. If the threat speed is s^, and 

the course is uniformly distributed over (0°,360°), the 

outer limit of the new planar region containing the threat 

after time t has elapsed, given that it was originally 

located in the p, * 100% confidence ellipse with semi-axes A 

and B, can be represented by an ellipse with semi-major axis 

A+s . t T and semi-minor axis B+s.t_ (Figure 7). Thus the 
1 Li 1 Li 
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FIGURE 7. TIME LATE PLANAR REGION 
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region over which the density will be evaluated is still 
the intersection of an ellipse with an angular wedge, 
recalling however that the density is no longer normal. 
Further, the new elliptical region does not represent a 
x 100% confidence region of the time late density. 
Calculation of the likelihood that the threat lies along 
bearing 3 given that the threat lies in the p^ x 100% ellipse 
at time zero, that the true bearing lies in (3 ,3 n ), that 
the sensor lies in the p, x 100% circle and that the speed is 
s ^ proceeds as follows: 

1. Estimate the parameters of the original normal dis- 

2 2 

tribution: mean = (0,0), = A /(-2 ln(l-p^)), 

o y 2 = B 2 /(-2 In ( 1 - p, ) ) . 

2 . Determine the upper and lower bounds on the true 

bearing wedge, 3 + k c; and g-k a Q . 

3 p p 3 

3. Determine the time late planar region as the ellipse 
with semi-major axis equal to A + s . t T and semi -miner axis 

X 1j 

equal to B + s . t . 

4. Determine the position of the sensor relative to the 
ellipse center. 

5. Subdivide the bearing fan ( 3 - k.a _ , 3 + k O o ) into 

3 P 3 3 

2n subintervals. 

6. At the midpoint of each subinterval, determine the 
equation of the line through the ship position at that 
bearing 3^. 

7. Let the equation of the line of bearing 3^ be 
Y-c 

X = . Then, the plane perpendicular to the xy-plane which 
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contains this line intersects the time late threat density 
in a curve whose equation is 



g L (y) 



= (■ 



m 



2 2 2 2 
a Y + m a x 






2 , 2 2 
°Y + m a X 



o — / 2 2 2* 

/2tt y a y + m a x 

2 2 2 

1 r / 'Y m J X w 

2 1 2 2 2 ^ 

m a x a Y 



v 2 rr m a a v 

A x 



c a. 



2,22 

ay + m a 



> 2 ] 



0 



360 . (s . t T ) cos ( 9-ip) -2s . t T ycos ( d-\p) 

f exp [- i 1 L ■ 1 L 



2 2 

m a x 



(s . t ) 2 sin 2 ( q—i/j) - 2s . t T ysin (8-ip) - 

1 L 1 xj ^ t a 

T 2 } 1 360 



(4) 



found by making the substitution x = in the density 

function f (x,v;t ) and expanding. 

A / X Li 

For computational purposes, assume that the width W 
of the region of the time late threat ellipse cut out by the 
angular subinterval around will be nearly constant through 
the ellipse. Let the points of intersection of the line of 
bearing 6^ with the time late threat ellipse be (X^^Y,^), 

( Xj. . , Y - ) . Then approximate the time late density at speed 
s. over this region by the absolute value of the product of 
and the area under the curve g between Y and Y.^ . The 
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area under curve g L can be approximated as follows. Sub- 
divide the interval (Y^, ,Y, 7 ) into n_ segments of length h. 
Evaluate g,. at the midpoint of each segment, y . Note that 
the first term in parenthesis in (4) is a constant equal to 

m times the value of the density of a normal random variable 

2 2 2 

with mean zero and variance a +m a v evaluated at c. The 

i A 

second term in parenthesis is the density function of a normal 
random variable with mean 



and variance 



cc, 



2,22 

i Y t m q ^ 



2 2 2 
m a x a y 

T~~ 2 

c y t m o x 



2 



evaluated at y. The variable y also appears in the integral 

term in equation (4). Numerically evaluate this term of 

(4) with y = y,. Let g_(y_) be the result of this computation. 

U i. vT 

Then, evaluating the area under the curve g T between Y^ 
and Y^.., is equivalent to the calculation: 



m 



lA 



2 J 2 2 
Y +m J X 



J = 1 






2,22 

C Y +m C X 






2 2 2 
a y +m o x 



X° Y 



(Z J> Si'i'j* 



where $ ( • ) is the density function of a standard normal random 
variable, and 
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C J, 



= - 



I Q 



2 2 2 
Qy + m 



2 : 2 2 
a + m a 

) ( : — — ) 

m °X°Y 



An analagous situation arises if the substitution Y = mX + c 
is made for y in the time late density f v (x,y;t_). 

A f X Xj 

8. The value of the density over the sub interval con- 
taining 3^ is then weighted by the instantaneous probability 
that the bearing error is 3^ - 3: 



W, 



k2 



f X, Y (X L' Y l'* t L ) dy L 



f„ ( 



■kl 



9. Steps C.7 and C.8 are then repeated for each sub- 
interval of (3 T ,B TT ) and the results of each calculation are 
Li U 

summed. 

10. The uncertainty in sensor location is accounted for 
by repeating steps C.4 through C.9 with the assumption the 
sensor is located at each of the three points in Figure 5, 
multiplying by the probability the sensor lies in that 
interval and summing each result. 

11. The result of calculations in steps C.l through 
C.10 is the likelihood that the threat lies on bearing 3 
given the threat was originally located in the p^ x 100% 
confidence ellipse, B_ e (3 T /3 rT ), the sensor is located in 
the p ? x 100% circle and the threat speed is s^. The con- 
dition that the speed is s ; is removed by repeating the 
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calculations C.l through C.10 for each of the speeds s^ , 
i = 1 , . . . ,M multiplying the result by the probability that 
the speed equals s. ; and summing all M results. The final 
result is the likelihood that the threat lies on 8 given 
it was originally located in the p^ x 100% ellipse, 

£ (3^/Sy), and the sensor is in the p^ * 100% circle esti- 
mated as 



M 



l P[S 
i=l 



s i ] 




( 3 ) . 
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IV. CONCLUSIONS 



Possible applications of and extensions to the algorithm 
developed in Chapter III are discussed in this chapter. 

As indicated in Chapter I, the objective of this paper 
has been to develop a means of assessing the likelihood a 
threat whose position has previously been estimated lies at 
a given measured bearing from a local sensor. The procedures 
were developed with a view towards permitting the user to 
make comparisons among lines of bearing measured by different 
sensors or among conflicting bearing information generated 
by one sensor. The approach chosen has been to estimate the 
likelihood that a threat lies on bearing 3 given that the 
threat is located in an ellipse of specified confidence 

x 100%, that 6 is measured with p. x 100% accuracy, and that 
the sensor position is measured with p^ * 100% accuracy. The 
algorithm to calculate the likelihood in the cases where the 
probability distribution of the target can be assumed to be 
bivariate normal at the time of the bearing measurement has 
been implemented on the TI-59 calculator. The cases in which 
this program applies are the following: (1) when the time 

elapsed since generation of the threat error ellipse is small 
enough to justify using the original estimate of the ellipse; 
(2) when threat course and speed are known, in which case 
the ellipse center is translated from the original position 
according to the course, speed and elapsed time information; 
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(3) when the threat can be assumed to be moving about in a 
random manner over a significant region in such a way that 
the ellipse center remains unchanged, but the x and y vari- 
ances have increased. When none of these cases hold, but 
the course is assumed uniformly distributed over (0°,360°) 
and speed has a discrete distribution over a finite interval, 
the threat density at the time late t T is not a bivariate 
normal density. The algorithm applicable in this case has 
not been implemented, but is described in some detail in 
Chapter III and Appendix B. 

Once the appropriate computation has been completed for 
each of the bearings considered, the results can be used to 
weight the value of several bearings in refining the threat 
location estimate provided by the external sensor. Note that, 
although unlikely, f'(3) may correctly be greater than one. 
Comparisons using these likelihoods should be made only when 
the upper and lower bounds on the true bearing fan for each 
bearing are chosen at the same probability level p„. and the 
uncertainty areas for all sensors include the same probability 
level p^. Further, these likelihood levels should be selected 
so as not to exclude a significant portion of the appropriate 
density. If p, or p, are not the same in all cases to be 
compared, f'(3) must be divided by the applicable value of 
p ^ or p. for each measured bearing 3 to be considered. 

Having established the relative value of available 
bearing information, the user can allocate weapons or further 
search effort accordingly. However, the probabilities calculated 
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are strictly ordinal data and do not define a redistribution 
of target location probability based on additional informa- 
tion. Further, the threat ellipse does not contain the 
target with certainty. The power to predict the probability 
of success of a search or weapons allocation plan based on 
the priorities established by these procedures is limited 
by these constraints. In this area in particular further 
research would be useful. 

In situations such as that for which the procedures in 
this paper have been developed, where a track has not been 
developed on the target, introduction of unrelated bearing 
data to a target motion model could impact significantly on 
the reliability of future position predictions. If there is 
high confidence in the reliability of the estimate of the 
threat ellipse provided by the ocean surveillance sensor, 
the prioritization established herein could be used to pro- 
cess bearing data prior to input to a target motion model. 

Using a pre-established threshold, only those bearings which 
coincide with the threat with an acceptable level of like- 
lihood could be used to refine or update a track on the 
threat . 

Desirable enhancements to the algorithms include providing 
for the instances in which the interval of uncertainty of 
the target course is known to be less than (0°,360°). Further, 
if the circular region of radius R contains the sensor with 
significantly less than 86% confidence or the assumption of a 
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bivariate normal distribution of sensor location is unsatis- 
cactory, it is left to the user to modify the calculations 
accordingly. 

The utility of these algorithms would also be improved 
by implementation on a larger and faster system than the 
TI-59 calculator. 
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APPENDIX A. TI-59 PROGRAM VERBAL FLOW AND USER'S INSTRUCTIONS 



Part I 
Step Number 

000 - 029 

020 - 029 
030 - 052 

053 - 053 
059-063 

064 - 073 



074 - 089 



Verbal Flow 

Enter the confidence level p, for 
the threat ellipse. Calculate the 
value of k for the given p : 

k = v -2 In (1-p, ) . 

Enter the length of the semi-major 
axis, A. Calculate o Y = A/k. 

Enter the length of the semi-major 
axis, B. Calculate a v = B/k. 

1 2 2 

Calculate B/A, aa v , a v and a v . 

X i a i 

Enter orientation of semi-major 
axis , e . 

Enter bearing from sensor to center 
of threat ellipse, a. Calculate 
9 — a . 



Enter distance r from sensor 
position to center of threat 
ellipse. Determine rectangular 
coordinates of sensor position 
( U, V) from polar coordinates 
(-r , 9-a) . 



Store the constants 360, /2tt . 
Initialize register 35 to 0. 
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090 - 096 



097 - 104 



105 - 143 



144 - 146 



147 - 151 



Enter number of standard devia- 
tions desired for bearing fan, k . 

3 

Enter standard deviation of bearing 
error, a,. Calculate k c,. 

Enter the angular stepsize desired, 
A3, for incrementally stepping 
through (3 L ,3 0 ). Calculate ^A3 . 
Calculate the largest number n of 
increments of size A3 contained in 
k,a. - jA3 degrees. Initialize 
counter 00 to n+1. Save n+1 in 
register 20. Determine 

5 A3 = k.,c, - iiS - nA3 / 

the residual increment. Calculate 

¥ 46 • 

Initialize counter 01 to 2. 
Calculations will be made at the 
midpoint of each interval from 
3 + i-AS to 3 + k a , then at the 

^ P P 

midpoint of each interval from 
3 - yA3 to 3 - k . j , , and finally 
at 3. Counter 01 indicates whether 
calculations are complete on both 
sides of 3 . 

Enter bearing measured by sensor, 

3 . 



152 - 154 



Enter standard deviation of sensor 
position a. 
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155 Enter index of the sensor position 

to be used for this run. 

156 - 157 Coordinates of the sensor position 

are selected in accordance with 
run number entered above in 
Subroutine sin. 

158 - 238 Determine whether a bearing 3' 

parallel to either axis is included 
in the fan (3 L ,3 lJ ). If 
includes a bearing parallel to the 
y-axis, use program 1 for Part II. 
The appropriate program number is 
displayed in calculator display 
register. If (3^,3^) does not 
contain a bearing parallel to 
either axis, use program 1. 



239 - 


259 


Subroutine 


P/R. 


260 - 


340 


Subroutine 


sin. 
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Part II, Programs 
Step Number 


1 and 2 

Verbal Flow 


000 - 003 


Initiate 3 ' . 


004 - 007 


If the last angular increment on 
this side of S has been considered, 
go to step 513 . Otherwise 
continue . 


008 - 012 


Decrement counter for angular 
increment. If counter = 0 , go to 
021. Otherwise continue. 


013 - 020 


Remove flag to indicate this is 
not the last angular increment. 
Recall A3/ the input angular step- 
size. Go to 029. 


021 - 028 


Set flag to indicate this is the 
last angular increment on this 
side. Add one-half the input 
angular stepsize A3 and one-half 
the residual stepsize 5 A3: 

^£3 + j5AB. 


029 - 030 


Increment 3 ' by the appropriate 
stepsize . 


031 - 045 


Convert 3' to an angle between 0° 
and 360° , 


046 - 051 


Calculate 9-3'. Print 6-3'. 
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052 



052 - 064 


If | 0 — 0 ' [ = 90° or = 270°, go to 236 


065 - 074 


If le-B'l = 0° or = 180°, go to 075 
Otherwise go to 133. 


075 - 086 


If the absolute value of the y 
coordinate of sensor position is 
greater than the length of the 
semi-minor axis, 3/ go to 004. In 
this case 3 1 does not intersect the 
error ellipse. Otherwise continue. 


087 - 088 


Set flag 2 to indicate that the 
bearing 6 ' is parallel to the 
x-axis . 


089 - 111 


Calculate the coordinates of the 
points of intersection of 3' with 
the threat ellipse. The y-coordi- 
nates are equal to V, the 
y-coordinate of sensor position, 
x-coordinates are determined in 
Subroutine v . The points of 
intersection are symmetric about 
the y-axis. Thus, X^ = -X . 

Store the smaller x value in 
register 27, the larger in register 
28. 


112 - 132 


Save locations of X,. , , X ^ ' 

a v , and a v . Go to 291. 

X i 


133 - 138 


Remove flags 2 and 3 to indicate 
that 3 * is not parallel to either 
the x- or y- axes . 
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139 



145 



146 - 154 



155 - 176 



177 - 235 



236 - 247 



Calculate the slope of s': 
m = tan ( 0 -3 ' ) . 

Calculate the y-intercept of 3 ' : 
c = V - mU. 

2 2 2 2 

If c > A'm" + B , go to 004. In 
this case, 6' does not intersect 
the error ellipse. Otherwise 
continue . 



Calculate and X, , the x- 

coordinates of the points of inter- 
section of 3 ' with the ellipse : 



-me + | + B" - c 2 

' kl = AV + B 2 



B r 2 2 '2 2 

-me - ~ y Am + B -c 

C k2 = “ 72 2 ” ^2 

Am +3 



Calculate the y-coordinates of the 
points of intersection of 3' with 
the ellipse, Y kl , 



Y ki = mX ki + c 



Y k2 “ mX k2 + C ' 



Go to 299 . 



If juj, the absolute value of the 
x-coordinate of sensor position, 
is greater than the length of the 
semi-major axis, A, go to 004. 3' 
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does not intersect the error 
ellipse in this case. Otherwise 
continue. 



248 - 249 



250 - 272 



273 - 298 



299 - 328 



Set flag 3 to indicate that 6 ' is 
parallel to the y-axis. 

Calculate the coordinates of the 

points of intersection of 3 ' with 

the error ellipse. The x- 

coordinates are equal to U, the 

x-coordinate of sensor position. 

The y-coordinates are determined 

in Subroutine y x . The points of 

intersection are symmetric about 

the x-axis. Thus, Y. , = -Y, 

kl k2 

Store the smaller y value in 
register 29, the larger in 30. 



Save the locations of Y^ , 
X kl ' °Y' a X* 



Determine the width of the strip 
around 3 ' , W: 



If flag 0 set, 

W = (d.+d 2 ) tan(j6AB), 

Otherwise 

W = (d^+d^) tan ^-AB) . 



dj and d. are the distances from 
(U,V) , the sensor position, to the 
intersection points ( X , 7 ) and 



( X k2 , Y k 2 ) » respectively . 



Distances 
are calculated in Subroutine log. 
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Multiply W by A3 if flag 0 not 
set. Otherwise, multiply by 6A3. 
Save result in register 37. 



329 - 334 



If 3' is parallel to either the 
x- or the y-axis, go to 468. 
Otherwise continue. 



When using program 1 ; 



335 - 380 



381 - 391 



392 - 401 



Express the y-coordinates of the 
intersection points as standard 
normal random variables, Z. and 



Z . 

l 



ca. 



‘ki 2/2 2 

J Y + * 3 X 

”, a X°Y 



T Jo" ~ + m a 



2 2 
X 



i = 1,2. 



Sort the values Z, in descending 
order. Let Z^’ be the larger 
value. Z ; ' is the smaller. 

Calculate 0(Z2*) - $(Z^'), the 
probability that a standard normal 
random variate lies between Z, ' 
and Z 2 ' . 
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402 - 424 



The result 



When using 
335 - 372 



373 - 383 



Multiply the results of steps 
392-401 by: 



/ F 2 2 

1/ a + m a 



*(■ 



z/cr „ 2 + m a 



•)x I m I x w. 



2 2 

X 



where 4> is the standard normal 
density function. 



of the calculations in steps 335-424 is 
Y, 



k2 

W / f 
*kl 



(Y — — y) dy 

X,Y k m ,y/ y 



program 2 : 



Express the x-coordinates of the 
intersection points as standard 
normal random variables Z, and Z . : 



Z . 

i 



me a. 



X k2 + 



J 2 ~ 2 2 

1/ ° v + m a 



a X a Y 



, / 2 ~ , 2 2 
1 / a x + m a Y 



i = 1,2. 



Sort the values of Z in descending 
order. Let Z 2 * = larger value. 

Z, 1 = smaller. 
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384 



393 



394 - 413 



Calculate 4> ( Z 2 ' ) - $(z 1 ') f the 
probability that a standard normal 
random variate lies between Z^' 
and Z 2 1 • 



Multiply the results of steps 
384-393 by: 



/ 



2,22 
°Y + m a x 






2 , 2 2 
°Y + m a x 



-) x W, 



where <p is the standard normal 
density function. 



The result of the calculations in steps 335-413 is 



W 




f (x,mx+c) dx 

A , X 



Regardless of which program is in use: 



425 - 465 



Multiply results of previous 
calculation by 



1 




where 




3 ’ 

6 '-360° 



3' < 180° 
S' > 180° 



and <j> is the standard normal density. 



The result of this calculation is 

f g (6;6') lw / f XrY (x,y) d Y 
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466 



467 



Go to 508 



468 - 507 



508-512 



513 - 515 



Calculate the area under the curve 
formed by the intersection of g ' 
with f y (x,y) . When g' parallels 
the x-axis this calculation 
becomes : 



(”-) I$(Z_) - $(Z. ) ] 

Oy 2 1 



where Z, 



X k2 / a x and S 1 X kl /a X 



V is the y-coordinate of the sensor 

position . 

When g' is parallel to the y-axis, 
the calculation is: 



*<!-) [HZ,) - -S(Z.) ] 
a X 3 X 2 1 



where Z. 



‘k2 //a Y 3 ld Z 1 



Y kl /3 Y 



U is the x-coordinate of the sensor 
position . 

Go to 422 to complete calculation 
of f ( 3 ; 3 1 ) | W /f x ^ y (x,y) | . 

(Go to 410 in Program 2) . 



Accumulate the probability at each 
angular interval. Display result. 

If flag 1 is set go to 558. 
Otherwise continue. 



516 - 520 If flag 0 is not set, that is if 

the angular increment just con- 
sidered was not the last on this 
side of 3, go to 004 and continue 
calculation on same side of 3- 
Otherwise continue. 
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521 



524 



525 - 542 



543 - 567 



558 - 566 



Decrement the counter 01. If 
counter is now 0 go to 543. In 
this case, the probabilities have 
been calculated on both sides of 3. 
The calculation at 6 remains to be 
done. Otherwise continue. 

Remove flag 0. Multiply A3, 4 ab, 

1 

and j6A3 by -1. Reinitialize 
counter 00 to n+1. Go to 000 to 
begin calculation on second side 
of 3 . 

Set flag 01 indicating that calcu- 
lations on both sides of 3 have 
been completed. The next iteration 
will do the calculation at 3' = S. 
Remove flag 00. Initialize 3 ' to 0 
Recall 3. Go to 029. 

Remove flag 01. Display the 
accumulated likelihood. STOP 



The result of this calculation is approximately f_ 1 (3) with 

sensor at (U.,V.) where i = run number entered in Part I. 

: 3 J 



567 - 587 



Subroutine y x calculates the points 
of intersection when 3 ' is parallel 
to either the x or y axis: 



x ki = - A i 1 - <l > 2 

X k2 = A lA "" S 
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Parallel to y-axis: 



589 - 610 



y ki ■ ' B tI 1 - <I> 2 

Y k2 - B A - (|> 2 • 

Subroutine log calculates the 
distance from sensor position to 
point of intersection of 8 ' with 
the ellipse: 

d i ■ !/Vx ki ) :; + <v-Y ki ) 2 
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part III 
Step Number 


Verbal Flow 


000 - 012 


Enter result of run 1. 
Multiply by 1/3. 


013 - 021 


Enter result of run 2. 
Multiply by 1/3. 


022 - 029 


Enter result of run 3. 
Multiply by 1/3. 


030 - 032 


Display likelihood. 


033 - 041 


If P 2 is the same for all bearings 
to be compared, enter 1. Go to 
042. Otherwise, enter p for this 
bearing. Divide likelihood by p-,. 
Display result. 


042 - 050 


If p^ is the same for all bearings 
to be compared, enter 1 and STOP. 
Otherwise, enter p., for this bearing 
Divide likelihood by p - . Display 
result. STOP. 
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USER'S INSTRUCTIONS 



The program to determine the likelihood that the threat 
lies along bearing ft given the threat is in the confidence 
ellipse, ftp e and the sensor is within the p^ x 100% 

confidence circle is in three parts. All parts require the 
use of a printer and the use of the Applied Statistics 
Library Module. Prior to running, the calculator must be 
repartitioned : 

1. Enter 4 

2. Press 2nd OP 17 

Part I 

1. Read sides 1 and 2 of Part I 

2. Read side 4 of Part II either program 1 or program 

2. Since program 2 is used more often, unless it is known 
that the bearing fan contains a bearing parallel to 

the major axis, recommend using program 2 of Part II. 

3. Press RST 

4. Enter p, ,, the confidence level of the threat ellipse. 
Press A 

5. Enter A, length of the semi-major axis. If data 
provided is length of the entire major axis, divide by 2 
before entering. 

Press B. 

6. Enter B, length of the semi-minor axis. If data 
provided is length of entire minor axis, divide by 2 before 
entering. Press C. 
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7. Enter 0 , the bearing of major axis from North. 

0° <9 < 180°. 9 is entered in degrees. Press D. 

8 . Enter a, the bearing of the threat ellipse center 
from the estimated sensor position. 0° £ a <_ 360°, in 
degrees. Press E. 

9. Enter r, the range from the estimated sensor posi- 
tion to the center of the threat ellipse. Press 2nd A'. 

10. Enter k , the number of standard deviations desired 

6 

in one direction from 3 . The program will construct a fan 
of equal size on the other side of 3 . Press 2nd B'. 

11. Enter a , the standard deviation of bearing error. 

p 

Press 2nd C 1 . 

12. Enter A3, the desired angular stepsize. Press 2nd 

D’ . 

13. Enter 3 , the measured bearing. Press 2nd E'. 

14. Enter a, the standard deviation of sensor position. 
Press R/S. 

15. Enter the number of this run, 1, 2 or 3. When run 
number = 1 , sensor is located at its estimated position 

(u ,v 0 ) . When run number = 2, the location will be 

O w 

(u o + .97/003(9-3-90) ,v + .97osin(0-B-9O) ) . When run number = 3, 
the location will be (Uq— . 97acos ( 0 - 3 - 90 ) ,v Q -. 97asin (© — 6~ 90) ) . 
Press R/S. If program number displayed matches that of the 
Part II side 4 read in, continue to 16. Otherwise, press 
2nd CMS, RST. Read in side 4 of the Part II program which 
corresponds to the number displayed. Repeat 3 through 15. 
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16. Press 4 2nd WRITE. Rerecord side 4 of the Part II 
program read in. This enables data entered in Part I to 
be transferred to Part II. 

Part II 

1. Read sides 1, 2, 3 and 4 of Part II Program 1 or 2 
as selected by the Part I program. 

2. Press RST 

3. Press R/S. Values of g and wff (g;g )f (g ) 

* 3 3 true 

will print alternately. Final result also prints out at end 
of calculation. 

4. Record final result: Likelihood threat lies along 

bearing g given threat is in p- ( x 100% ellipse , true bearing 

is in p- x 100% fan and sensor is at location used in this run. 
Parts I and II must be completed 3 times (Run numbers 1, 2 
and 3) before proceeding to Part III if uncertainty in sensor 
position is being considered. 

Part III 

1. Read side 1 of Part III. 

2. Enter result of run 1. Press A. 

3. Enter result of run 2. Press B. 

4. Enter result of run 3. Press C. 

5. Enter p 0 , confidence level of (g ,3 ), if necessary. 

Z Li U 

Otherwise enter 1. Press D. 

6. Enter p., confidence level of sensor position, if 
necessary. Otherwise, enter 1. Press E. 
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APPENDIX B. TIME LATE VERBAL FLOW 



VERBAL FLOW 



1 . 


Enter 


P 1 ‘ 


Determine the value of k: 


k = /-2 In 1 


2 . 


Enter 


A. 


2 2 

Determine 0 = (A/k) 

A 




3. 


Enter 


B. 


2 2 

Determine = (B/k) 




4. 


Enter 


9, 


orientation of major axis 




5. 


Enter 


a , 


bearing from sensor to ellipse center 


6 . 


Enter 


range r from sensor to ellipse 


center 


7. 


Calculate 


(u ,v ) the coordinates of 
0 0 


mean sensor 



position. u Q = -r cos(Q-a), v Q = -r sin(Q-a) 

8 . Enter a , the standard deviation of bearing error 

3 

9. Enter k , the number of standard deviations to be 
included in the bearing fan on each side of 3 . 

10. Enter desired angular stepsize, AS 

11. Calculate the number of iterations of size AS required 

on each side of g: I = [ (k a - 4 a B ) /AS ] where [•] means 

S3 ^ 

the greatest integer less than or equal to the value within. 

12. In general (k.a,. -i-A3)/A3 is not an integer. Determine 

3 p ^ 

the size of the fractional increment: 

6 AS = { (k_ - - ~-AS)/AS - I) AS • 

13. Enter 3, the measured bearing. Let PROB2=Q, PROB1=0 

14. Enter the time late t 

1j 

15. For each of the discrete threat speeds to be con- 
sidered, repeat steps 16 to 43. Then go to 44. 

16. Enter the target speed s. Let PROB = 0 
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17. Expand the outer limit of the threat ellipse: 



Let A 



A + st. 



L 



B + st 



L 



18. Let 8' = 6 

19. For each of the I increments of size A3 repeat 
steps 20 to 38. 

20 . Let 3' = 3 ' + AS 

21 . Calculate (9-3') . If j 6 — B ' j = 0° or 180°, go to 27. 
If | 0-6' | = 90° or 270°, go to 30 

22. Calculate m = tan (9-3')/ the slope of the line of 
bearing . 

23. Calculate c = mU+V, where (U,V) is the sensor loca- 
tion for this iteration. c is the y-intercept of the line 
of bearing. 

2 2 2 2 

24 . Calculate m A +B -c . If this quantity is less 
than zero, the line of bearing s' does not intersect the 
threat ellipse. Go to step 20 and process next increment of 
size AS/ if any remain. If all I intervals of size AS on 
this side of 3 have been considered go to step 40. 

25. Calculate the points of intersection of the line 
of bearing with the expanded ellipse: 



X 



1 




X, 



2 
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mX^ + c 



Y 1 = mX ± + c, Y 7 = 

26. Go to 32 

27. If | V | > B', the bearing 6' parallel to the x-axis 
does not intersect the threat ellipse. Go to step 20 and 
process next increment of size A3 if any remain. If all I 
intervals if size A3 on this side of 3 have been considered, 
go to 40. 

28. Calculate the points of intersection 

y i - y 2 - V 

X ± = -B ' jjl - (V 2 /B ' 2 ) 

X 2 = B'/l - (V 2 /B' 2 ) 

29. Go to 32 

30. If | U | > A', the bearing 3' parallel to the y-axis 
does not intersect the threat ellipse. Go to step 20 and 
process the next increment of size A3 if any remain. If 
all I intervals of size A3 on this side of 3 have been con- 
sidered, go to 40. 

31. Calculate the points of intersection: 



Y 1 = -A' t/ 1 - (U 2 /A' 2 ) 

= A' /l - (U 2 /A' 2 ) 
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32. Calculate the median width of the strip of the 
ellipse defined by the angular subinterval under considera- 
tion. If the subinterval is of size A 6 : 

di + d 2 ^ 

W = -=~r — — ( 2 tan ±A3) 



If the subinterval is of size 5 A 



a l u 2 1 

W = — — -(2 tan joAS) 



In these expressions d^ is the distance from the sensor to 
the ith point of intersection, i = 1, 2: d = 1 / (U-X^) ' + (V-Y^)~ 

33. If | 6 - 3 ’ 1 = 0° or 180°, go to 36. If j 0-3 ' j = 90° 
or 270°, go to 35. 

34. Evaluate the target density from Y, ' to Y„ } 1 along 
the line y = mx + c: 

Let Y 1 ' = min (Y^Y,) 

Y 2 ' = max(Y 1 ,Y 2 ) 

<)3 (z) = density function of a standard normal 

random variable evaluated at z 



K = (m//j Y ~ + m 2 a x 2 ) 4 (c//a y 2 + m 2 c x 2 ) 



Subdivide the interval ( Y ^ , Y „ ' ) into n segments of length 
h. At the midpoint of each segement, y ^ , compute: 
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2 2 

+ m a x / ma a ) ) 



360 

/ exp[- 




2 v 2 2 



0 



+ 




when the integral must be numerically evaluated. 

Sum the results over all n segments. 

Go to 37. 

35. Evaluate the target density from 1 to 1 along 
the line x = U: 

Let Y x ' = min ( Y ,Y , ) 



Subdivide the interval (Y into n segments of length 



Y 2 ' = max(Y 1 ,Y 0 ) 



.j; (z) = density function of standard normal random 



variable 



K 
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where the integral must be numerically evaluated. 

Sum the results over all n segments. Go to 37. 

36. Evaluate the target density from 1 to X 0 ' along 
the line y = V: 

Let X^ ' = min^-^X^) 

X, 1 = max(X^,X ) 

( z ) = density of standard normal random 
variable 



Subdivide the interval (X ' ,X., ') into n segments of length 
h. At the midpoint of each segment, x ^ , compute: 



. x. 360 
hK — ( <j> (-^-) ) 

J X a X 0 



(stj^) cos (d-\p) - 2stLx. (cos (e-ip) 

exp [- j ( — - — 2 3 

C X 

(st^) 2 sin 2 (e-’p) - 2st i ysin(e-^) 



360 



where the integral must be numerically evaluated. Sum 
results over all n segments. 

37. Multiply the value of the target density just 
computed by W. 

38. Multiply this result by 



1 



<p( 




where <P is the standard normal density function. 
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39. Let PROB = PROB+ (the results of the calculations 
in steps 20 through 38 for each if the I subintervals of 
size 13) • 

40. If the fractional interval of size 6A3 on this 
side has been considered go to 41. Otherwise, let 

3' = j6A3 + jA3- Repeat calculations 21 through 38 once. 
Let PROB = PROB+ (the result calculated at this step) . 

41. If the computations on both sides of 3 have been 
computed, go to 42. Otherwise repeat the computations 
from 18 to 40 on the other side of 3 by letting A3 = -A8/ 

5 A 3 = -5A3- 

42. Let 3' =3. Repeat steps 21 through 38 once. 

Let PROB = PROB+ (the results of this calculation) . 

43. Let j = the number of discrete threat speeds to 
be considered. Let PROB1 = PROB1 + P[S=s]PROB. Go to 15. 

44. Repeat steps 15 through 43 once with the sensor 
located at each of three points: 



(1) 


(U,V) = 
density; 


(u ,v ) , the mean of the sensor 
o o 


(2) 


(U,V) = 


(u q + . 9 7ocos (0-3-90) , 
v q + .97asin(e-3-90)); 


(3) 


(U,V) = 


(u - .97 cos (3-3-90) , 
v q - . 97 sin (9-3-90) ) ; 



where a is the standard deviation of the sensor density. 
The value of PROB1 calculated at each iteration will be 
weighted by the approximate probability that the sensor is 
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located in the region of the sensor error circle repre- 
sented by the applicable value (U,V) . If p is of the 
order .86 or greater, multiply by 1/3. 

45. The value of PR0B2 calculated after completion 
of step 44 is the relative probability that the threat 

lies along bearing 3 at time t , given the threat was in 

Li 

p. x 100% ellipse at time t , Q. is in (3 T /S TT ) and the 
1 o t L U 

sensor is in a p^ x 100% confidence circle. 
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CALCULATOR PROGRAM 
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35 
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42 
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04 
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43 
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06 
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43 
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10 
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43 
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06 


06 
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42 


STD 
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07 


07 
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043 


42 
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045 


43 
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